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Abstract 

Let M and N be compact smooth oriented Ricmannian rt-manifolds 
without boundary embedded in K n+1 . Several problems about min- 
imal distortion bending and morphing of M to N are posed. Cost 
functionals that measure distortion due to stretching or bending pro- 
duced by a diffeomorphism h : M — > N are defined, and new results 
on the existence of minima of these cost functionals are presented. In 
addition, the definition of a morph between two manifolds M and N 
is given, and the theory of minimal distortion morphing of compact 
manifolds is reviewed. 



1 Introduction 

Two diffeomorphic compact embedded hypersurfaces admit infinitely many 
diffeomorphisms between them, which we view as prescriptions for bending 
one hypersurface into the other. We ask which diffeomorphic bendings have 
minimal distortion with respect to some natural bending energy function- 
als. More precisely, let M and N be diffeomorphic compact and connected 
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smooth oriented n-manifolds without boundary embedded in R n . The 
manifolds M and ./V inherit Riemannian metrics qm and and correspond- 
ing volume forms from the usual metric and orientation of R n+1 . Although 
we will use this structure here, only the existence of the metrics on M and iV 
is essential. We pose the problem of bending M into ./V via a diffeomorphism 
h : M — > N so that the distortion produced by h is minimal with respect 
to some functional that measures bending or stretching (cf. the problem of 
optimal development of surfaces [16j [7] ) . 

The problem of minimal distortion bending of manifolds may be con- 
sidered a special case of the problem of minimal morphing. A morph is 
defined to be a transformation between two shapes through a set of in- 
termediate shapes. A minimal morph is such a transformation that mini- 
mizes distortion. There are important applications of minimal morphing in 
manufacturing 0116], computer graphics OQj)], movie making [XT] , and 
mesh construction [HE]. We will formulate and solve a problem about the 
existence of minimal morphs with respect to stretching for n-dimensional 
manifolds. 



2 Distortion Minimal Bending 

Let (M, <7Af) and (N,gw) be smooth compact Riemannian n-manifolds with- 
out boundary with all the additional properties stated in the introduction. 
By Vol(M) we denote the volume of M. Also, let Diff (M, N) denote the set 
of all diffeomorphisms between M and N. 

Our first functional measures distortion due to stretching. For a point 
p G M, we define the distortion due to stretching at p as the infinitesimal 
relative change of volume. More precisely, let {Ak}^ =1 C M be a sequence 
of open neighborhoods of the point p that shrink to the point as k — > oo. 
For example, one can choose A% = B(p, x) D M, where B(p,R) is the open 
ball of radius R in W n+l centered at p G M C R n+1 . 

Definition 2.1. The distortion due to stretching produced by a diffeomor- 
phism h G Diff (M, N) at a point p G M is defined to be 



fh(A k ) 


iV N \ 













i{p) = lim ^ - , fc = J(h)(p) - 1, 

fe^oo 

where J(h) is the Jacobian of h with respect to the (Riemannian) volume 
forms ojm and uijsf. The functional <3?i : Diff (M, N) — > M+ is defined by 

^(h) = [ (\j(h)\ -iYlum. (1) 
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The following results are proved in [3]. 
Lemma 2.2. A diffeomorphism h G Diff (M, N) is a critical point of 3>i if 



and only if J(h){m) = y°)|^ for all m 



G M. 



Theorem 2.3 (Existence of minimizers for <I>i). If [M,qm) and (N,g^) are 
diffeomorphic compact connected oriented Riemannian n-manifolds without 
boundary, then there exists a minimizer of the functional 3>i over the class 
Diff(M, N) and the minimum value of is 

$min = (Vol(M)-Vol(AQ) 2 
1 Vol(M) 

The functional $1 is invariant with respect to compositions with volume 
preserving maps: <3?i(/i o k) = 3>i(/i) provided that k G Diff(M) is volume 
preserving (it has the Jacobian J(k) = 1). Therefore, the minimizer of $1 
is not unique. 

For a vector bundle V over M, we denote by T(V) the space of all sections 
of V; T r (V) denotes the space of all C r sections of V. Let T^ 2 \M) denote 
the vector bundle of covariant order- two tensors over M (see pQ). In order to 
measure distortion with respect to bending, we introduce the strain tensor 
field S : Diff(M, N) -> r(T(°' 2 )(M)) by S(h) = h*g N - g M , where h*g N is 
the pull-back of the metric g^ by h. 

Definition 2.4. The deformation energy functional $2 : Diff(M, N) — ► R+ 
is given by 

|2, 



$ 2 {h) = / \\h*g N -9m\\ mm, 
Jm 

where the fiber norm || • || on the bundle T^ ' 2 \M) is induced by the fiber 
metric G := g* M ® g* M (see [12]). 

Let r°°(TM) denote the space of C°° sections of the tangent bundle of 
M. Fix a diffeomorphism h G Diff(M, N). In order to derive the Euler- 
Lagrange equation, we consider all variations of the diffeomorphism h of the 
form ho(j) t , where <pt is the flow of a vector field Y G T°° (TM) . Because the 
tangent space Diff (M, N) can be identified with the space of all sections 
T°°(h~ 1 TN) of the pull back bundle h _1 TN over M, every smooth variation 
of the diffeomorphism h can be represented in the form hofa (see pj for a 
more detailed description). 

The diffeomorphism h G Diff (M, N) is a critical point of <&2 if 



d_ 
~dl 



$2(Wi)| t= o = £>$ 2 (/t)^ = 2 / G(h*g N -g M ,L Y h*g N ) = (2) 

Jm 
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for all Y € r°°(rM), where h*Y is the push forward of Y by h, and Ly is 
the Lie derivative in the direction of Y (see [1] ) . 

Let V be the Riemannian connection on M generated by the Riemannian 
metric qm with Christoffel symbols T^- (see |10|). The connection V with 
Christoffel symbols f ^- is the Riemannian connection of the metric h*g^ on 
M. 

Definition 2.5. Define B(h) = (h* g^ — gu)* '■ That is, for each p e M, the 
tensor B(h)(p) of type (2,0) is defined as the tensor (h*gjy — <?m)(p) with 
its indices raised. 

Definition 2.6. Define the bilinear form A(h) to be 

A{h)(X,Y)=V x Y -V X Y (3) 
for X,Y e r°°(TM) (see [12] ). 

It is easy to prove that A(h) is a tensor field of type (1,2) on M with 
components 

W hp kp fcp- 

Lemma 2.7. T/ie /irsi variation of the functional $2 ^ i^e direction Y € 
r°°(TM) is given by 

D$ 2 (h){KY) = -4 I g M (d\vB(h)+A(h) : B(h),Y)u M , (4) 

J M 

where B(h) : A(h) is the contraction of the tensor fields B(h) and A(h) (see 
Moreover, h is a critical point of the functional $2 if and only if 

dwB(h) + A(h) : B(h) = 0. (5) 

The Euler-Lagrange equation for the functional $2 is the system of nonlinear 
partial differential equations ©. 

Let hji : M n+1 — > M n+1 be the radial map given by hji(x) = Rx for some 
number R > and for all x £ M n+1 . It is easy to check that if N = RM is 
a rescaling of the manifold M, then the map h = hn o f satisfies the Euler 
Lagrange equation ([5]), whenever / G Diff(M) is an isometry on M. 

The following results on minimizing $2 in the one-dimensional case are 
proved in [2]. 

Proposition 2.8. (i) Suppose that M and N are smooth simple closed 
curves in M? with arc lengths L(M) and L(N) and base points p G M and 
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q € N; 7 and£ are the arc length parametrizations of M and N with^{Q) = p 
and £(0) = q that induce positive orientations; and, the functions v and w 
are defined by v(t) = L(N)/L(M)t and w(t) = -L(N)/L(M)t + L(N) for 
all t e [0, L(M)}. If L(N) > L(M), then the functional <J> 2 has exactly two 
minimizers in the admissible set 

A={he Diff (M, N) : h(p) = q} : 

the orientation preserving minimizer 

h\ = £ o v o 

and the orientation reversing minimizer 

h-z = £ o w o 

(where we consider^ as a function defined on \0,L(M)) so that^/~ l {p) =0). 
Moreover, the minimum value of the functional <3?2 is 

_ {L{N? - L{Mff 
$2 " L{Mf • (6) 



Proposition 2.9. Assume the notation of the proposition \2.8l 

(i) If L(N) < L{M), then the functional $2 has no minimum in the admis- 
sible set 

Q = {h £ C°°(M, N) : h is orientation preserving and h(p) = q}. 

(ii) If jjj^ < 4| 3 then the functional $2 has no minimum in the admissible 
set A = {h € Diff (M, N) : h{p) = q}. 

The main ingredients for a proof of (i) in proposition 12.91 are simply 
illustrated. The curve M is wrapped around the curve N without stretching 
and the excess is removed. This wrapping function can be expressed in the 
form h = £cuo7 -1 (in the notation of proposition [2T8|) , where u : [0, L(M)] — > 
R is a discontinuous piecewise linear function. The function h is not smooth; 
but, it is possible to approximate it by a minimizing sequence {h k }^L l C Q 
whose deformation energies $>2(h ) converge to &2{h) = 0. On the other 
hand, ^(Z) > for all / 6 Q. The proof of (ii) uses the second variation 
of <J? 2 - 

By proposition 12.91 we see that even in the one-dimensional case the 
functional $2 exhibits nontrivial behavior: the minimum does not always 
exist, and the existence depends on properties of the curves M and N. 
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The general problem of the existence of minimizers for the functional <E>2 
is open. On the other hand, we have solved the problem for the case where M 
and iV are Riemann spheres or compact Riemann surfaces of genus greater 
than one. Let H(M,N) = {h G Diff(M, N) : h is a holomorphic map}. 

Theorem 2.10. (i) Let Iir : M 3 — > M 3 be the radial map given by hji(p) = 
Rp for some number R > 0. // M = S 2 C M 3 and N = h R (M), then 
h := / o Iir is a global minimum of the functional §2, restricted to the 
admissible set 7i(M,N), whenever f is an isometry of N . 

(ii) Let M and N be compact Riemann surfaces. If H(M, N) is not 
empty and the genus of M is at least two, then there exists a minimizer of 
the functional <E>2 in Ti(M, N). 

The general problem of minimization of the functional $2 seems to be 
very difficult because the admissible set is an infinite-dimensional manifold 
Diff(M, N) whose structure is not completely understood. In theorem 12. 101 
the admissible set is a finite-dimensional homogeneous space in case M and 
N are two-spheres and a finite group in case M is a compact Riemannian 
manifold of genus greater than one. A natural idea is to reformulate the 
problem of minimal distortion bending in such a way that the admissible set 
is a linear space. 

Fix a diffeomorphism / G Diff(M, N). Every diffeomorphism h : M — > 
N can be represented in the form h = f o (ft, where eft € Diff(M). For 
simplicity of notation, let g<i := f*gN an d gi ■= gu- To measure the de- 
formation produced by h = f o <p £ Diff(M, N), we use the strain tensor 
field S(4>) = (f)*g2 — 9i- In other words, the problem reduces to minimiza- 
tion of the deformation energy produced by some class of diffeomorphisms 



The tangent bundle TM is equipped with the Riemannian metric g\. 
Let W k,2 (TM) be the (k, 2)-Sobolev space of sections of the tangent bundle 
TM (see [13J ) . We choose the number of generalized derivatives k G N 
large enough so that the Sobolev space W k,2 (TM) is embedded into the 
space r 2 (TM) of all C 2 sections of TM and some additional estimates hold. 
Consider the space H = L 2 ([0, 1]; W k ' 2 {TM)) of time dependent vector 
fields v : M x [0, 1] — > T(TM). The space H is a Hilbert space equipped 
with the norm 



Every vector field v G H generates a diffeomorphism on M in the fol- 



0: (M )fl i) -> (M,g 2 ) in Diff(M). 
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lowing sense. The nonautonomous ordinary differential equation 



-dt =V{q > t) 



(7) 



is solved (on the compact manifold M) by an evolution operator r] v (t;s,p) 
that satisfies the Chapman-Kolmogorov conditions (see [3]) and is such that 
rj v (s;s,p) = p for every p £ M. The function (f) v : M — > M given by 
4> v {p) = rf(l\ 0,p) is called the time-one map of the evolution operator if; 
it is a diffeomorphism on the manifold M. 

We define the distortion energy functional E : H — > M + to be 



where Ilj is the second fundamental form on M associated with i = 1,2 
(see |10j). This functional incorporates strain (which is intrinsic to the 
manifold M) and bending (which is extrinsic). 

Theorem 2.11 (Existence of minimizers for E). There exists a minimum 
of the functional E in the space H. 

The proof uses the direct method of the calculus of variations as well as 
convergence properties of the evolution operators if (t;s,x) generated by 
weakly convergent sequences of time dependent vector fields in H. 

3 Distortion Minimal Morphing 

Definition 3.1. Let M and ./V be compact connected oriented n-dimensional 
smooth manifolds without boundary embedded in M n+1 . A C 1 function 
F : [0, 1] x M — > M ra+1 is a morph from M to N if the following conditions 
hold: 

(i) p i— > F(t,p) is a diffeomorphism onto its image for each t £ I = [0,1]; 

(ii) the image M t = F(t, M) is an n-dimensional manifold possessing all 
the properties of M and iV mentioned above; 

(hi) p i— > F(0,p) is a diffeomorphism of M; 

(iv) the image of the map p F(l,p) is N. 



We denote the set of all C 2 morphs between the manifolds M and N by 
M(M,N). 




gi\\ 2 UM + 



JM 



[ ||OTn 2 -iii|| 2 u, M , (8) 
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For simplicity, we will only consider morphs F such that p \—> F(0,p) 
is the identity map. We assume that each manifold M = F(t, M) (with 
M° = M and M 1 = N) is equipped with the volume form u% = Lfl, 
where Q, = dx\ A dx2 A ... A dx n+ \ is the standard volume form on R n+1 
and rjt ■ M t — > M ra+1 is the outer unit normal vector field on M* with 
respect to the usual metric on Also, as a convenient notation, we use 

/* = F(t, •) : M — > MK 

Definition 3.2. The functional : Diff (M S ,M*) -> M + is defined by 
formula {!]), were M and iV are replaced with M s and M* respectively. A 
morph F is distortion pairwise minimal (or, for brevity, pairwise minimal) 
if / s '* = /* o (/ s )~ x : M s — > M* minimizes the functional for every 
s,t £ [0,1]. We denote the set of all C 2 distortion pairwise minimal morphs 
between manifolds M and N by VM{M, N). 

Using lemma [2721 it is easy to derive a necessary and sufficient condition 
for pairwise minimality. 

Proposition 3.3. Let M = M° and N = M 1 be n-dimensional manifolds 
as in definition \3.1\ equipped with the ( respective ) volume forms loq and lo\ . 
A morph F between M and N is distortion pairwise minimal if and only if 



for all t € [0, 1] and m £ M , where J(f ) is the Jacobian of f with respect 
to the volume forms ujq and uit- 

The following proposition states the existence of pairwise minimal morphs. 
It can be proved by rescaling morphs between M and N, which are not nec- 
essarily pairwise minimal, to conform to property (jUJ). Moser's theorem on 
volume forms (see [9]) plays a crucial role in the proof. 

Proposition 3.4. Let M and N be n-dimensional manifolds as in defini- 
tion \S.l\ If M and N are connected by a C 2 morph, then there is a distortion 
pairwise minimal morph between them. 

Having the preliminary study of pairwise minimal morphs at hand, we 
define minimal morphs. 



Definition 3.5. The infinitesimal distortion of a C 2 morph F from M to 
N at t G [0, 1] is 



■/(/*) (m) _ 1 



(9) 



Vol(M*) Vol(M) 



e F (t) = lim 




M 



s 



where E s,t = 3>i'*(/ s '*) is the distortion energy of the transition map f s,t . 
The total distortion functional 9 defined on such morphs is given by 



i f i 



The following proposition implies that it suffices to minimize the func- 
tional 9 over the class VAi(M, N) of all pairwise minimal morphs instead 
of the class M(M, N) of all C 2 morphs. 

Proposition 3.6. (i) The following inequality holds: 

inf 9(G) < inf 9(P). (11) 

GGVM(M,N) P&M(M,N) 

(ii) If there exists a minimum F of the total distortion functional 9 over 
the class VM(M, N), then F minimizes the functional 9 over the 
class Ai(M, N) as well; in fact, 

9(F) = min 9(G) = min 9(P). (12) 

G&VM(M,N) P&M{M,N) 

Using proposition 13,31 it is easy to recast the functional 9 into a simpler 
form: 

Lemma 3.7. The total distortion of a C 2 pairwise minimal morph F from 
M to N is 

The latter form of the functional 9 and proposition 13.61 allow us to solve 
the problem of minimization of the total distortion functional 9 over the 
class of all C 2 morphs A4(M, N). In order to solve the problem, we minimize 
the auxiliary functional 

m = J o 1 jdt (14) 

over the admissible set 

Q= {4>e C X ([Q, 1];K+) : 0(0) = Vol(M),0(l) = Vol(JV)}. 

The following theorem — our main result on distortion minimal morphing— 
is proved using proposition 13.61 and lemma [3771 

Theorem 3.8. Let M and N be two n-dimensional manifolds satisfying the 
assumptions of definition \3.1[ If M and N are connected by a C 2 morph, 
then they are connected by a minimal morph. The minimal value of 9 is 

min 9(F) = A(JVo\(N) - JVol(M)) 2 . (15) 

f<eM(M,N) v ; 
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